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We propose a solid-state experiment to study the pro- 
cess of continuous quantum measurement of a qubit state. 
The experiment would verify that an individual qubit stays 
coherent during the process of measurement (in contrast to 
the gradual decoherence of the ensemble-averaged density ma- 
trix) thus confirming the possibility of the qubit purification 
by continuous measurement. The experiment can be realized 
using quantum dots, single-electron transistors, or SQUIDs. 



The impressive advantages promised by quantum com- 
puting have revived the interest to the fundamen- 
tal quantum effects in simple objects: two-level systems, 
which in this context are nowadays called qubits. In this 
paper we address the problem of continuous measurement 
of a qubit state having in mind a solid-state realization 
of the setup. 

Among the numerous proposals of quantum comput- 
ers, the solid-state realizations (see, e.g. Refs. (|-||]) 
look more promising because of better controllability of 
qubit parameters and inter-qubit couplings. However, 
the qubit measurement in this case is not as straightfor- 
ward as in typical optical experiments where the single 
photon just "clicks" the detector. The reason is finite 
(and typically weak) coupling with a solid-state detector 
and finite intrinsic noise of the detector. As a result, the 
measurement cannot be done instantaneously, and so the 
collapse postulate of the "orthodox" quantum mechanics 

cannot be applied directly. Instead, the quantum mea- 
surement should be considered as a continuous process. 

There are two main theoretical approaches to the con- 
tinuous quantum measurements. One approach (which 
dominates in solid-state physics and so can be called 
"conventional" ) is based on the theory of interaction with 
dissipative environment 0,^. Taking trace over the nu- 
merous degrees of freedom of the detector, it is possi- 
ble to obtain the gradual evolution of the density ma- 
trix of the measured system from the pure initial state 
to the incoherent statistical mixture, thus describing the 
measurement process. Since the procedure implies the 
averaging over the ensemble, the final equations of this 
formalism are deterministic and can be derived from the 
Schrodinger equation alone, without any notion of the 
state collapse. 

The other approach (see, e.g., Refs. |s|-[l^]) is closer 
to the collapse viewpoint and describes the stochastic 
evolution of an individual quantum system due to con- 



tinuous measurement. This evolution obviously depends 
on a particular measurement result and is usually called 
selective or conditional quantum evolution. Depending 
on the details of the studied measurement setup and ap- 
plied formalism, different authors |P-p^ discuss quan- 
tum trajectories, quantum state diffusion, stochastic evo- 
lution of the wavefunction, quantum jumps, stochastic 
Schrodinger equation, complex Hamiltonian, method of 
restricted path integral, Bayesian formalism, etc. The 
theory of selective quantum evolution was only recently 
introduced into the context of solid-state mesoscopics 
In particular, it was shown that the continu- 
ous measurement of an individual qubit does not lead 
to gradual decoherence (in contrast to the conventional 
result for the ensemble), instead, the measurement can 
lead to gradual purification of the qubit density matrix. 

Since the concept is still considered controversial, the 
experimental check is quite important. In this paper 
we propose an experiment which can be realized using 
three possible setups available for present-day technol- 
ogy: double-quantum-dot qubit measured by quantum 
point contact JlGf , qubit based on single-Cooper-pair box 
measured by single-electron transistor |17(], or SQUID- 



based qubit measured by another SQUID |18| , |19[ . 

Let us start with reviewing the result of the conven- 
tional formalism for the continuous measurement (Fig. 1^) 
of a qubit state (see, e.g. recent publications ]l^ , p0| -p5 1 ) . 
For the qubit characterized by the standard Hamiltonian 
^QS — (£/2)(c|ci — c\c2) + H{c\c2 + eld) in the basis 
defined by coupling with the detector, the evolution of 
qubit density matrix p is given by equations 



pii = -p22 = -2i7Impi2, 

pi2 = iepi2 + iH{pii - P22) - rpi2. 



(1) 
(2) 



where the continuous measurement is described by the 
dephasing rate F, which was calculated for different mod- 
els in Refs. [|6|,|oH||]. 

These equations do not depend on the detector output 
because they represent the result of ensemble averaging, 
including the averaging over the measurement result. To 
study the evolution of an individual qubit let us denote 
the noisy detector signal as I{t) (assuming current for 
definiteness). Two "localized" qubit states 1 and 2 cor- 
respond to average detector currents Ji and I2 which by 
assumption do not differ much, A J = Ji — /2 ^ -^0 = 
(Ji + /2)/2. (This assumption of "weakly responding" 
detector allows us to use the linear response the- 
ory and also Markov approximation if the processes in 
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FIG. 1. Schematic of a qubit continuously measured by a 
detector with output signal I{t). 



the detector are much faster than the qubit evolution.) 
Intrinsic noise of the detector signal is characterized by 
the spectral density S which is frequency-independent in 
the range of interest. The noise determines the typical 
measurement time t ^ S/ (A/)^ necessary to distinguish 
between states 1 and 2, and thus defines the timescale of 
the selective evolution of the qubit density matrix p{t). 

Within the Bayesian formalism the selective evo- 
lution is described by equations 

pn = -2mmpi2 + (2 pnP22 [I{t) ~ h] , (3) 
Pl2 = «£Pi2 + iH{pii - P22) 

-(A//5)(pii -p22)[/(<) -/o]pi2 -7P12, (4) 

where the dephasing 7 = F — (A/)^/4S' > is now due 
to the contribution from "pure environment" only. In 
particular, 7 = if the qubit is measured by symmetric 
quantum point contact, since in this case F = (AI)'^ /AS 
(see Refs. (20|j2|Jl|,||l). We will call such detector an 
ideal detector, 77 = 1, where 77 = 1 — 7/F is the ide- 
ality factor. In contrast, the single-electron transistor 
p6[ in the operation point far outside the Coulomb 
blockade range is a significantly nonideal detector [p4| , 
77 <C 1; however, rj becomes comparable to unity when 
the current is mostly due to cotunncling processes [p7| . 
The SQUID is an ideal detector when its sensitivity is 
quantum-hmited [|8|j2||. 

Eqs. allow us to calculate the evolution of qubit 

density matrix p if the detector output I{t) is known from 
a particular experiment. To simulate the measurement 
we can use the replacement H] 



det. A 




qubit 




det. B 



I{t)-Io = AI{pn-p22)/2 + at), 



(5) 



where the random process ^{t) has zero average and 
"white" spectral density = S. One can check that 
averaging of Eqs. over all possible measurement 

results [i.e. over random contribution ^(t)] reduces them 
to Eqs. (|l|)-(||). Notice that the stochastic equations are 
written in Stratonovich form which preserves the usual 
calculus rules, while averaging is more straightforward in 
Ito form 

As follows from Eqs. (|^)-(^), if a qubit with initially 
pure state, |/9i2(0)P = /5ii(0)/322(0), is measured by an 
ideal detector, then its density matrix p{t) stays pure 
during the measurement process. Even if initial state is 
a statistical mixture, p{t) is gradually purified during the 
measurement 
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The predictions of the Bayesian formalism can be 
checked experimentally, however, it is not quite simple 
at the present-day level of solid-state technology. The 



FIG. 2. Schematic of two-detector correlation experiment 
using Cooper-pair box and two single-electron transistors. 



direct experiment was discussed in Ref. The idea 

was to perform the measurement by almost ideal detec- 
tor during some finite time r, record the detector output 
I(t), use Eqs. (||)-(|^) to calculate p{t) and then check 
the calculated value. This check can be done by chang- 
ing qubit parameters e and if in a way to ensure pu = 1 
at some specified moment of time, that can be measured 
by the detector switched on again. Since for coherent 
evolution the qubit can be placed with 100% certainty in 
the state 1 only if the wavefunction is known precisely, 
such check (repeated many times) verifies that p(r) is 
pure and coincides with the calculated value. 

Unfortunately, this experiment would require very fast 
recording of I(t). Since the expected coherence time is on 
the order of 10-100 ns at most (see, e.g. |l^), the band- 
width of the detector signal coming out of the cryostat 
should be at least 1 GHz, that is very difficult experimen- 
tally. Another proposed experiment |2^,0 is to measure 
the spectral density of the quantum coherent oscillations 
and check the predicted maximal peak-to-pedestal ra- 
tio of 4. Such an experiment may be easier to realize 
(because the basic spectral analysis can be done on-chip 
inside the cryostat), however, it would not prove unam- 
biguously the Bayesian formalism, since an alternative 
interpretation of the result is possible |^ . 

Here we propose an experiment which is even easier to 
realize, and which can test the Bayesian formalism (|^)- 
(^). The main idea is to use two detectors {A and B) 
connected to the same qubit (Fig. ^). The detectors are 
switched on for short periods of time by two shifted in 
time voltage pulses (one for each detector) with durations 
TA and tb, supplied from the outside. The output signal 
from the detector A is the total charge Qa = Jq^ I a {t) dt 
passed during the measurement period. Similarly, the 
output from the detector B is Qb = J^'^'^" Isit) dt, 
where r is the time shift between pulses. If the mea- 
surement by the detector A changes the qubit density 
matrix, it will affect the result of measurement B. Re- 
peating the experiment many times (with the same initial 
qubit state) we can obtain the probability distribution 
P{Qa,Qb\t) of different outcomes, which contains the 
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information about the effect of the quantum measure- 
ment on the qubit density matrix. In comparison with 
previous suggestions, the advantage of this correlation 
experiment is that the wide signal bandwidth is required 
only for input pulses (that is relatively simple) while the 
outputs are essentially low frequency signals. The experi- 
ment can be called "Bell- type" because of some similarity 
with the famous proposal of Ref. |^ . 

Fig. 1^ shows the realization of the experiment using 
single-electron transistors (two small tunnel junctions in 
series [^) as detectors. Qubit is realized by the Cooper- 
pair box |3^,|l^ so that the electric charge of the central 
island can be in coherent combination of two discrete 
charge states. Another similar setup is two quantum 
point contacts measuring the charge state of a double- 
quantum-dot qubit. One more setup is the 3-SQUID ex- 
periment in which the qubit is realized by one SQUID 
while two other SQUIDs are in the detecting regime. For 
definiteness we will consider the realization of Fig. ||. 

The conventional formalism ([l])-(H) does not give any 
explicit predictions for the resulting probability distribu- 
tion P{Qa,Qb\t)- However, it implies the absence of 
correlations between p{t) and I{t), so for example the 
average result of the second measurement Qb{Qa,t) = 
J QbP{Qa,Qb\t) dQs should not depend on Qa- The 
Bayesian formalism (^)-(^ makes the different predic- 
tion: does depend on Qa- 

For simplicity let us assume symmetric qubit, £ = 0, 
which is initially in the ground state, pn — P22 = 
P12 — 0.5, and also assume relatively strong coupling 
between the qubit and detectors, (AIa)^ /HSa 3> 1, 
{AIb)'^ / HSb ^ 1 (subscripts A and B correspond to 
two detectors), so that we can neglect the qubit evo- 
lution due to finite H during the measurement periods 
TA and tb, which are assumed to be on the order of 
Sa.b/{AIa.b)^. Then from Eqs. (|)-(|) if follows that 
the first measurement only "partially" localizes the qubit 
state and after obtaining the result Qa from the first 
measurement the qubit density matrix is 



{Qa - taI2a)'^ - {Qa - taIia)" 



2/511 (ta) - 1 tanh 

2bATA 

Pi2{ta) = [pii{ta) P22{ta)]^^^ ejip{--fATA) , 

where Eq. (||) is actually the classical Bayes for- 
mula which was used in Ref. jl^ to derive the for- 
malism (|3|)-(^). [The probability to get Qa has 
the distribution P{Qa) ~ {pi + ?'2)/2 where pi ~ 
{ttSata)'^^'^ e:iq){-{QA - taIia)'^ / Sata)-] The qubit 
performs the free evolution during the time t — ta be- 
tween measurements (here we neglect ta <?C t) and 
the average result of the second measurement Qg — 
tb{I2B + Pii{t)A.Ib) depends on Qa in the following 
way (Fig. 3a): 

^ 1 . {Qa-taI2A? -{Qa-taIiaY 
= 2 2S^A 
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FIG. 3. The normalized average result (5s of the second 
measurement for several selected results Qa of the first mea- 
surement, as a function of the time r between measurements. 
Panels (a)-(b) are for strong coupling and panels (c)-(d) for 
moderate coupling between the qubit and detectors (other 
parameters are the same). The calculations are done by 
Bayesian formalism while the conventional formalism does not 
predict any nontrivial dependence. 
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X cos ( rir — arcsin cxp(— 7/t/2) , (8) 

where Sb = {Qb ^ tbIob)/tbA.Ib, 7/ is the dephas- 
ing with both detectors switched off, and fl = {4H^ — 
7^/4)^/^ is the frequency of quantum oscillations (under- 
damped case is assumed). Notice that Sb changes sign 
together with the sign of Qa — taIqa, while the phase of 
oscillations is a piece-constant function of Qa- 

The dependence becomes quite different if the 7r/2 
pulse is applied to the qubit immediately after the first 
measurement, that multiplies Pi2{ta) given by Eq. (|^) by 
the imaginary unit. In this case (Fig. 3b) 

(6) 5b= Asm{flT + arcsin z / A) exp(— 7/r/2) , 

(7) A^[{z' + y'-yz^f/2H)/{l-^j/16H')]'/\ (9) 

where z = Pii{ta) — 1/2 and y = Impi2(r4 -t- 0) = 
Repi2{TA — 0) are given by Eqs. (||)-(@). This expression 
considerably simplifies for weak dephasing, jata <C 1 
and 7/ <ti H , when 



1 Tf 
(5b = - sin[fiT -I- arcsin(2pii(rA) - 1)] exp( 



(10) 



In contrast to Eq. (^ now the phase of oscillations 
(5b (r) depends on the result Qa of the first measurement, 
while the amplitude is maximal possible and independent 
of Qa- This fact is very important since it proves that 
after the first measurement (by an ideal detector) the 



3 



control 
stage 



feedback 
signal 



desired p(t) 



comparator 



31 



qubit ^ — ) detector > eqs. (3)-(4) 



I(t) 



P(t) 



environment 



FIG. 4. Schematic of the continuous qubit purification 
using the quantum feedback loop. 



qubit remains in the pure state for any result Qa- This 
state depends on Qa and is not one of the locaUzed states 
as somebody could naively expect. [Notice that Eq. (||) 
can in principle be interpreted in terms of such "classical" 
localization, as indicated by its independence on r]A-] 

In a realistic experimental situation the assumption of 
strong coupling with detectors may be inapplicable. In 
this case the full probability distribution P{Qa, Qb\t) as 
well as the dependence Qb{Qa,t) can be calculated nu- 
merically using Eqs. The results of these calcu- 
lations for (AIa^/HSa = iAlB)^/HSB = 1 are shown 
in Figs. 3c and 3d. Weak coupling as well as the non- 
ideality of the detectors decrease the correlation between 
the results of two measurements, however, for moderate 
values of the coupling and nonideality the correlation is 
still significant. 

The successful experimental demonstration of the cor- 
relation and quantitative agreement with the results of 
the Bayesian formalism would prove the validity of this 
formalism and therefore prove its other predictions. Be- 
sides the clarification of the relation between the mea- 
surement result and qubit evolution, the important for 
practice prediction is the gradual qubit purification due 
to continuous measurement which can be useful for a 
quantum computer. 

All quantum algorithms require the supply of "fresh" 
qubits with well-defined initial states. This supply is not 
a trivial problem since the qubit left alone for some time 
deteriorates due to interaction with environment. The 
usual idea is to use the ground state which should be 
eventually reached and does not deteriorate. However, to 
speed up the qubit initialization we need to increase the 
coupling with environment that should be avoided. The 
other possible idea is to perform the projective measure- 
ment after which the state becomes well-defined. How- 
ever, in the realistic case the coupling with the detector 
is finite that makes projective measurements impossible. 
Here we propose a different way: to tune qubit continu- 
ously in order to overcome the dephasing due to environ- 
ment and so keep qubit "fresh" . 

The idea of such state purification is shown in Fig. ^. 
The qubit is continuously measured by weakly coupled 
detector, and the detector signal is plugged into Eqs. (^)- 
which allow us to monitor the evolution (in particular, 
due to interaction with environment) of qubit density ma- 



trix Pij{t). This evolution is compared with the desired 
evolution and the difference is used to generate the feed- 
back signal which controls the qubit parameters H and 
e in order to reduce the difference with the desired qubit 
state. We have performed the Monte-Carlo simulation of 
the qubit purification by feedback loop (in the regime of 
well-pronounced quantum oscillations) and found strong 
suppression of the qubit dephasing due to environment 
in the case when the dephasing rate 7 is comparable or 
weaker than the "measurement rate" (A/)^/4S'. 

In conclusion, we have proposed the Bell-type experi- 
ment which can test the Bayesian formalism predictions 
for the evolution of an individual qubit due to continu- 
ous quantum measurement. The next (and much more 
difficult) step is the experimental realization of the qubit 
purification using quantum feedback loop. 

The author thanks L. P. Rokhinson, D. V. Averin, M. 
H. Devoret, C. M. Marcus, and K. K. Likharev for useful 
discussions. 
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